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We generalize the treatment of inflationary perturbations to deal with the non-Markovian colored 
noise emerging from any realistic approach to stochastic inflation. We provide a calculation of 
the power-spectrum of the gauge- invariant comoving curvature perturbation to first order in the 
slow-roll parameters. Properly accounting for the constraint that our local patch of the Universe 
is homogeneous on scales just above the present Hubble radius, we find a blue tilt of the power- 
spectrum on the largest observable scales, in agreement with the WMAP data which show an 
unexpected suppression of the low multipoles of the CMB anisotropy. Our explanation of the 
anomalous behaviour of the lower multipoles of the CMB anisotropies does not invoke any ad-hoc 
introduction of new physical ingredients in the theory. 
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I. INTRODUCTION 

During the last two decades the inflationary paradigm 
has become an almost universally accepted scenario to 
explain the observed large scale flatness and homogene- 
ity of the Universe [HQ- At the same time, it pro- 
vides an efficient mechanism (the stretching of quantum 
fluctuations up to cosmological scales) to generate those 
small curvature perturbations (of the order of 10 -5 ) that 
are believed to produce the observed Cosmic Microwave 
Background (CMB) temperature anisotropies and to pro- 
vide the seeds for the formation of the cosmic structures 
(galaxies and clusters) we observe today. 

Stochastic inflation represents one possible ap- 
proach to the inflationary paradigm, specifically con- 
ceived to describe the transition from quantum fluctua- 
tions to classical density perturbations |4| . The basic idea 
is to introduce a cutoff in Fourier space through a suit- 
able time-dependent window function that filters out the 
modes whose frequency is lower than the comoving hori- 
zon size. The inflaton field is thus split in two differently 
behaving parts: the short-wavelength part has an in- 
trinsically quantum nature, while the coarse-grained one 
(collecting the remaining super-horizon modes) is treated 
as classical. A Langevin-like equation of motion for the 
long-wavelength part is obtained, where the sub-horizon 



'Electronic address: michele.liguori@pd.infn.it 

^Electronic address: matarrese@pd.infn.it 

t Electronic address: marcello.musso@pv.infn.it 

§ Electronic address: antonio.riotto@pd.infn.it 



modes enter now as a classical stochastic noise term that 
perturbs the dynamics of the coarse-grained field. 

In most works on this subject |20l | the noise is taken 
to be a white noise, whose two-point correlation function 
is by definition proportional to a Dirac delta function in 
time. A remarkable feature of stochastic processes involv- 
ing a white noise, known as Markov property, is that their 
conditional probability does not depend on times before 
the constraint. For this particular kind of noise, stan- 
dard techniques of statistical physics p| allow to write 
an evolution equation (the Fokker-Planck equation) for 
the probability distribution of the inflaton perturbations. 

A key issue in the solution of the Langevin or Fokker- 
Planck equation is the choice of the initial conditions 
for the perturbations. Many authors agree that it 
should be consistent to assume the spatial homogeneity of 
our observable local patch of the Universe, and therefore 
the vanishing of all fluctuations right before the moment 
it crosses the horizon size, about 60 e-folds before the 
end of inflation, since at that time only fluctuations on 
larger scales could have grown up. Therefore, all points 
inside the present Hubble radius (at that time contained 
in the same coarse-graining domain) must have the same 
local value of the scalar field, although this value can be 
different from the one assumed in other regions of the 
Universe. Even if it is generally assumed that inflation 
started well before the last 60 e-folds, for the white-noise 
case the evolution of fluctuations is completely insensitive 
to what happened before that epoch and the constraint 
really becomes a new initial condition. 

However, a white noise arises only if the window func- 
tion is chosen to be a step-function in Fourier space, 
therefore if coarse-graining is achieved with the introduc- 
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tion of a sharp momentum cutoff. Even though it must 
be regarded as the mathematical limit of a smoother fil- 
tering procedure, it has been pointed out @ that the 
step function is a somewhat pathological choice, because 
it does not satisfy some fundamental properties charac- 
terizing "well-behaved" smooth window functions. Con- 
versely, all these window functions always produce a non- 
divergent colored noise with the same asymptotical be- 
haviour. Unfortunately, the statistical treatment of a col- 
ored noise is a rather more difficult matter , but in gen- 
eral it is expected that because of their non-Markovianity 
fluctuations will keep some memory of the evolution be- 
fore the constraint. 

In a previous paper we studied the dynamics of 
the fluctuations of the coarse-grained inflaton field in 
a pure de Sitter space-time and found evidence for a 
blue tilt in the power-spectrum on the largest observ- 
able scales as a consequence of the non-Markovian dy- 
namics near the constraint. This is due to the fact that 
the increased noise correlation time (with respect to the 
white-noise case) acts as a sort of "inertia" against the 
growth of the perturbations after the constraint, thereby 
resulting in a suppression of the power-spectrum on the 
scales that crossed the horizon in the following few Hub- 
ble times. 

This is an interesting feature, since the CMB 
anisotropy measurements made by WMAP [llj give some 
evidence for a suppression of the low multipoles, con- 
firming earlier analogous results found by COBE [l2^ . 
Although the statistical significance of such a suppres- 
sion is not large [l3[ . many authors have recently tried 
to explain this spectral feature either invoking astrophys- 
ical effects |T3 | or introducing some new physical input 
in the mechanism that generates the perturbations [lfij . 

In this paper, we address the question of the low mul- 
tipoles suppression from the stochastic inflation point 
of view, suggesting that this might be simply a conse- 
quence of the colored noise. Compared to a white noise, a 
smooth choice of the window will in fact slightly suppress 
the contribution to the noise given by the field modes 
whose frequency is immediately higher than the cutoff 
scale o~(aH) (while enhancing lower frequencies) where a 
is a parameter smaller than unity, which is introduced to 
parametrize to level of arbitrariness - besides the shape 
of the window function - in the size of the coarse-graining 
domain. Right after the time r* at which we set the ho- 
mogeneity constraint on our comoving patch of Universe, 
fluctuations with k < oa^H* will grow less than in the 
white-noise case before freezing out, and if a is not too 
small this suppression can be effective also on observable 
scales. However, there seems to be no physical motiva- 
tions to state that the coarse-graining domain has to be 
very large compared to the horizon size. This assump- 
tion indeed appears in early stochastic inflation works, 
where a is taken to be much smaller than unity, but 
this is done essentially for practical reasons, in order to 
ease up the calculations and to cast away from observable 
scales the possible effects introduced by the white-noise 



approximation. This choice is justified a posteriori with 
the argument that the noise correlation function in con- 
figuration space is independent of a up to second order. 
Nevertheless, a a dependence still remains (even in the 
Markovian case) in the power-spectrum. 

The plan of the paper is as follows: in Section ITT1 we 
briefly review the formalism of stochastic inflation, ex- 
tending it to the colored-noise case [13 ■ In SectionlTTIlwc 
exactly calculate the colored-noise correlation function 
for a free massive scalar field in a de Sitter space-time. 
In Section HVI we evaluate the effects of the constraint at 
t* computing the conditional probability distribution of 
the same scalar field. In Section[V]we extend our results 
to the perturbations of the inflaton field (including in our 
treatment also metric fluctuations) and we compute the 
power-spectrum Vn of the gauge-invariant curvature per- 
turbation, taking as the only approximation a first-order 
slow-roll expansion of the background inflaton field. Fi- 
nally, our conclusions are drawn in Section fVTl 



II. EFFECTIVE STOCHASTIC ACTION 

We begin our considerations by taking a free scalar 
field with mass m in a homogeneous and isotropic back- 
ground space-time, whose metric has the pure Robertson- 
Walker form 



ds 2 = dt 2 - a 2 (t)dx 2 = a 2 (r)(rfr 2 - dx 2 ), 



(11.1) 



where a(t) is the scale factor responsible for the expansion 
of the Universe and r is the conformal time, defined in 
such a way that dr = dt/a(t). 

The ordinary action for a scalar field in a curved space- 
time with this background metric reads 



S[4>] = d A xa 3 (t) 



(V0) 5 



m <p 



(11.2) 



and the equation of motion for the normal modes 4>k (as- 
suming a spatially flat Universe these are simply obtained 
by expanding the field on the 3-dimcnsional Fourier basis 



„zk-x 



/(2tt) 3 / 2 ) is 



0k + 3ff0 k + —0k + m 2 0k = 0, (II.3) 

where dots denote ordinary time derivatives. 

It is convenient to define a new field variable x> whose 
normal modes are related to the Fourier modes 0k by 



Xk = 0ka(r); 



(II.4) 



in conformal time, indicating with a prime the conformal- 
time derivative d T = adt, the equation of motion for the 
new field becomes 



Xk = 0, 



(II.5) 
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and can be obtained by variation of the conformal action 



S\x] = ~ J dix \ X 



□ + mV 

a 



X, 



(11.6) 



where the metric determinant does not appear and □ = 
q2 _ y2 - 1S ordinary flat space-time d'Alcmbcrt oper- 
ator. 

In the de Sitter case, the Hubble parameter H = a/a 
measuring the expansion rate is strictly constant in time, 
and the scale factor evolves as 



a{t) 



Mi 



(II.7) 



while the conformal time, ranging from — oo to 0, be- 
comes 



a(t)H' 



(II. 



The evolution equation for the Xk's can then be recast in 
the form 



Xk 



1 



1 



t' \ 4 
where the parameter v is defined as 



Xk = 0, 



9 m 2 



4 H 2 ~ 2 6 ' 



(II.9) 



(11.10) 



The generic solution to this equation, expressed in 
terms of Bessel functions of the first and second kind, 
is 



fll.ll) 



c iy /\r\Mk\T\)+C2y/\T\Y v {k\T\); 



requiring each Xk to match the plain wave solution 
e~ lkT /\/2k for k ^> aH, when wavelengths are too short 
to feel any space-time curvature effects, produces the 
standard Bunch-Davies solution 



Xk(r) 



Hi 1] {k\r\), 



where 



HW(x) = J v (x)+iY v (x) 



(11.12) 



(11.13) 



it the Hankel function of the first kind. For the normal 
modes of the original scalar field <fi, this gives 



Mr) = ^H\r\^H^(k\r\), 
that in the massless case (v = |) becomes 



Mr) = H 



kr 



—ikr 



(11.14) 



(11.15) 



Each mode is said to leave the horizon when the corre- 
sponding physical wavelength afc -1 encompasses, due to 



the growth of the scale factor a(t), the size of the cosmo- 
logical particle horizon H . At the horizon crossing we 
therefore have \kr\ = 1. 

It is possible to obtain the sub-horizon part of the field 
selecting the short-wavelength modes, through a suitable 
time-dependent high-pass filter in Fourier space. We do 
this by means of a window function W a (kr) such that 
W a (kr) = for k\r\ < a and W a (kr) = 1 for k\r\ > a, 
where the parameter a (in early stochastic inflation works 
called e) introduces some level of arbitrariness in the size 
of the coarse-graining domain, corresponding to the "ef- 
fective horizon" a(aH). The fine-grained sub-horizon 
field </>> is therefore defined as 

W a (kr) 



dk- 



a k k (r)e 



-ik-x 



h.c]. 



(11.16) 



(2tt) 3 / 2 

while the coarse-grained super-horizon part cf> < = Lp is 
left for the moment unspecified. 

An effective equation of motion for the coarse-grained 
part can simply be obtained by substituting cj> = tp + 0> 
directly in the equation of motion (|II.3J| j| ; all the high- 
frequency part is then collected into a random noise 
field which acts as a classical source term for the long- 
wavelengths part. The quantum fluctuations on sub- 
horizon scales are averaged into a classical noise term 
£ (with a given probability distribution P[£\) perturbing 
the super-horizon dynamics: the effects of such a per- 
turbation is then a stochastic process to be studied with 
the ordinary methods of classical statistical physics. The 
quantum problem of computing the expectation value of 
the coarse-grained field <p is thus reduced to the classical 
problem of evaluating the mean of the solution to the 
stochastic evolution equation averaged over all possible 
noise configurations. 

However, things are a bit more complicated than this, 
since we are dealing with true expectation values between 
the same in state, and not simply with in-out transi- 
tion amplitudes (which are the usual objects in ordinary 
quantum field theory). The calculation of expectation 
values is a typical problem of out-of-equilibrium field 
theory [Tfil |. In this framework, an effective action for 
(tp) can be obtained via the so-called influence functional 
method 0, where the splitting of the field is operated 
in the action, and the short wavelengths are integrated 
out with a path-integral over all the configurations of the 
sub-horizon field (j) > . This method actually introduces 
some extra terms in the effective equation of motion jlDj . 
However, these effects are small and are neglected in the 
present work. 

Since we take the effect of sub-horizon fluctuations on 
large scales to be small, we can split the field into its 
statistical mean value tp (satisfying the classical equation 
of motion (|II.3|> ~) and a fluctuation <fy>[£], that by defini- 
tion vanishes when averaged over the £'s. The stochastic 
equation of motion for the super-horizon fluctuations can 
be reduced to 



5tp\ + 3HStp w — — - m )5tpk 



Ck 

n3- 



(11.17) 
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or equivalcntly 



k 



+ m 2 a 2 



(11.18) 



which are the usual equations obtained in stochastic infla- 
tion inserting the frequency split of the field directly into 
the equation of motion. This is a fully classical Langevin- 
likc equation, where the noise £ (resulting from the av- 
erage of the quantum fluctuations on sub-horizon scales) 
stochastically perturbs the super-horizon dynamics, act- 
ing as a source for the long-wavelengths field ip, which in 
turn can be treated classical stochastic variable. 

The noise £ is a Gaussian random field, whose config- 
urations are weighted by the functional probability dis- 
tribution 



P[Z) 



N exp 
iVexp 



d 4 xd 4 x'£{x)A- 1 (x,x')Z{x') 



(11.19) 



-\ Jdrdr'dkdk'U^A^ir, t')&(t') 
where Ar!,(T,r') is the functional inverse of 
A k ,k'(r, t' 



Re[/(fcT)/*(fcrQ] 
(5(k + k) — , (11.20) 



and 



(11.21) 



This probability distribution allows us to calculate the 
statistical mean value (■■■}§ of any ^-dependent quantity 
averaged over all the noise field configurations, defined as 



(...) s = Jv[t]...P[£\. (11.22) 

Then, by definition the mean (£(t)) s of the noise van- 
ishes at all times, while the two-point correlation func- 
tion, which completely characterizes the statistical prop- 
erties of the Gaussian noise field, is then by definition 



(&(r)&/(7 y )) s = A k)k ,(t > f) (11.23) 
In configuration space the correlation function reads 



s = 



dk 



3 ik(x— x ) 



'\ 1 



s '(2tt) 3 " 2fc 3 

Rc[/(fcr)/*(fcr')]. (11.24) 

The statistical behaviour of the noise thus critically 
depends on the shape of the filter. Choosing the usual 
window Wcr(fcr) = #(fc|r| — a) we obtain the standard 
white-noise two-point correlation function. For x = x' it 
reads 



Hi 

4tt 2 



(l + 0(a 2 ))S(t-t'); (11.25) 



which is highly divergent for t = t' and has a vanishing 
characteristic correlation time, while a smooth window 
yields a correlation function with no divergences and a fi- 
nite correlation time, therefore producing a colored noise. 
Namely, choosing 



W a {kr) = 1 - e"^ 3 



(11.26) 



the two-point correlation function for r = can be cal- 
culated, yielding 



mat')), 



i 



in 2 coslT (H(t-t')) 



0(a 2 ), (11.27) 



that asymptotically behaves like e~ 2H (* -t ). Moreover, 
it is possible to show that this asymptotic result is quite 
general for a wide class of smooth window functions Q ■ 



III. FLUCTUATIONS 

The particular solution of the evolution equation 
(|II.18(I for the fluctuations Sxk. sourced by the noise field 
£ can be expressed in terms of the general solutions 
Xi = ^k\r\J v {k\T\) and X 2 = y/ie\f\Y v (k\T\) of the ho- 
mogeneous equation i|II.9|) . This solution reads 



*Xk[f](r)= / d?g(kT,kf)&(f), 



(III.1) 



where 



^^) = ^i5'^^ (ni.2) 

and Ti n is the beginning of inflation, at which we set the 
initial condition Sxk(Ti n ) = 0. 

Keeping this assumption, we now want to introduce in 
our system the constraint that at a much later time t* 
(roughly about 60 e-folds before the end of inflation) we 
have no fluctuations in that part of Universe correspond- 
ing to the present observable sky. This is motivated by 
the fact that in our treatment all the points we observe 
today (over which we measure a substantial homogene- 
ity) were included at in the same coarse-grained do- 
main. It is thus consistent to assume t* the complete ho- 
mogeneity of the comoving patch of Universe we observe 
today, all fluctuations on smaller scales being generated 
later by the noise term. In this approach, we conserva- 
tively make no assumptions on the behaviour on larger 
unobservable scales. 

We are thus led to consider (for a given noise configu- 
ration) a different solution for the subsequent evolution 
of the fluctuations, obtained as in 
integration at 



MI. If) by starting the 
when a new (stochastic) initial condi- 
tion holds. In turn, i$Xk[£]( T *) is determined again from 
(|III.1|I with the usual vanishing initial condition at n n . 
However, as far as we are dealing with points inside the 
present observable Universe, we can skip the stochastic 
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initial conditions t* since their inverse Fourier transform 
is assumed to vanish. Therefore, in configuration space 
the subsequent evolution of the fluctuations will only con- 
tain noise modes integrated after r*. That is, for relevant 
x's we write 



function lTL20t. We find 



5ip(x, t) = 
and 



IWW~ 2 ~ J r f9(kT,kf)Uf) (ni.3) 



dk 



-J dfg(kr r , kf) £ k (f) 



(2tt)3/ 2 a* 

(III.4) 

where the first equation is only valid for scales inside our 
observed patch of the Universe. 

As expected, since the fluctuation <5y>k[£] is linear in £, 
at all times we have that 



(M£](t)) s = o, 



(III.5) 



while the two-point correlation function in xi and X2 = 
xi + r can be obtained integrating the noise correlation 



J 



C(r) S (^K]( Xl ,r)^K](x 2 ,r)), 



(27T) 



2k 3 



(III.6) 



where 
h(k) 



J T dfg(kr,kf)(W^ X k + 2<xL)- ( IIL7 ) 



In order to evaluate this integral, we first change the 
integration variable from f to x = kf and apply the 
relation 

J v (x)-^Y v (x) - ^J v (x)Y v (x) = — (III.8) 
ax ax irx 

in the denominator of g(kr,kf). Second, we integrate 
by parts the W'J term and and use again l|III.8|) . The 
integral can thus be written as 



h(k)=H 



-xW' a {x)[Y v {k\r\)Ju{x) 



J v {k\r\)Y v {x) H v {x) 



k\r,\ 



Hu{k\T\) 



k\r,\ d 

dx — W a (x) 
m dx 



(III.9) 



r 



containing only a boundary term and a trivial integra- 
tion. 

In the same way we can calculate the correlation func- 
tion evaluated at t», yielding 



C4r)^M£](xi^*)M£](x2,T*)>, 



dk 



\h{k)Y 
2fc 3 



(111.10) 



where hik) has the same form as h{k) but it refers to 
the time interval [t^t*]. However, the boundary term 
is now proportional to W^(/c|rj„|) and its contribution to 
the correlation function will be effective only on the scales 
that crossed the horizon at the beginning of inflation. 
Therefore, since in most realistic models r* S> i~; n , the 
boundary term in our treatment is completely negligible, 
and wc have 



h{k) = H 



7rfc 3 | 



-H v (k\r.\) 



(W a (k\r in \)-W a (k\n\)y (111.11) 



Moreover, we can also define the (generally non- 
vanishing) mixed correlation function 



M(r) = (ty[f](xi,r)ty[fl(x a ,T.)> s 
dk_ eikr R*[/i(fc)/2*(*0] 



(2tt) 3 



2fc 3 



(111.12) 



involving the scalar field perturbations evaluated at dif- 
ferent times. 



IV. PROBABILITY DISTRIBUTION 



We can construct the joint probability distribution 
P\6<pi, 5(f2, c] that the stochastic variables Jy[£](xi) and 
(5y[^](x2) satisfy the constraint c (i.e. they vanish) at 
time r* and assume the values dtpi and Stf2 at time r by 
taking the statistical average 
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P[5<p 1} 6<p a ,c] = ( [] Sfa -<M£](*^))*(Mf](*i>'-.)J ) (IV.l) 

\i=l,2 / s 

= y daida 2 d0 1 dP 2 c _ iy a . Sv> . J a<M£](*»r)-M2: /3«<5¥>K](x«,T„)p£]_ 

Since the field is linear in the noise £, the evaluation of the mean consists in a simple Gaussian functional 

integration yielding (repeated indices are summed) 

P{5<p 1} 5<p 2 ,c) = / ^ e -i:^,.-|^(r,) t . 3 f « e -jAO.( r «)ft-o (M (,«)A i (IV . 2) 



where rn = r22 = while ri2 = r 2 i = r, and thus 
the result now only involves 2x2 symmetric correlation 
matrices. 

According to Bayes theorem, the conditional probabil- 
ity P c (S<pi, S(p 2 ) is the joint probability P(S<p\, Sip 2 , 0, 0) 
normalized by the probability -P*(0, 0) of the constraint. 
The latter can be evaluated exactly in the same way, tak- 
ing the mean value over the noise configurations of the 
product of two (5-functions constraining the value of the 
fluctuations at t*. Following the same steps as before, 



the constraint probability reads 

P*(c)=/^e-i^(^ 1 (IV.3) 

and thus provides the correct normalization needed in or- 
der to evaluate in ljiV.2|) the Gaussian integration over /3i 
and (3 2 . With these results, we are now able to compute 
the conditional two-point correlation function 



P(Sip 1 ,Sip 2 ,c 
(<fy>(xi)(fy(x 2 )) c = / 5ip 1 Sip 2 - 



P*{c) 
J (2tt) 2 

M(0)M(r) C,(0)(M(0)-M(r)) 2 

C(r) " 2 a(o) + a(r) + c?(o)-c?(r) ■ (IV - 4) 



This result looks rather messy, but it considerably sim- 
plifies if we take the very reasonable limit t» 3> Ti n (as it 
is the case in most inflationary models). Actually, we can 
assume that I 2 {k) is given by i|III.9|) (with r* instead of 
r) neglecting the boundary term, since this would mod- 
ify the correlation function only on extremely large scales 
(those that crossed the horizon at the beginning of infla- 
tion). We then get 

~ (k\T*\y(w a (k\r m \) - PMfclr.i)), (iv.5) 

which is non- vanishing only for It^I" 1 < k < |t»| _1 , and 
since the spatial fluctuations are effective for k > r _1 > 
| t - * | 1 and are therefore damped by the cutoff introduced 
by Wo^/cIt*!) we roughly have 



which is divergent for |Tj„| — * oo and e — > 0. 

However, the dependence on Tj„ contained in I 2 (k) 
is rapidly saturated in the mixed correlation function: 
actually, it enters in M (r) only through the product 
I\(k)I 2 (k), and I\(k) (which is again obtained from 
(|III.9|I with instead of Tin) significantly differs from 
only for k > | -r^ | 1 . Therefore, when integrating over % 
the modes responsible for the r^-dependence (i.e. those 
such that k ~ |ri„| _1 ) are completely damped, M(r) is 
constant in the limit Ti n <C t» , and in this approximation 
we get 

(Mxi)ty(x 2 )) c ~ C(r). (IV.7) 

The power-spectrum Ps<p(k) of the fluctuations, de- 
fined so that 

<Mxi)M*a)>c = ^ / dk e* kr ^^, (IV.8) 
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becomes 



Air 2 



(IV.9) 



For relevant scales and small values of a we can 
gleet the boundary term in 



ne- 

l|III.9(l . In the small-tr limit 
we therefore recover the standard scale-invariant result 
V Sv {k) = H 2 /4tt 2 . 



V. CURVATURE PERTURBATIONS 

We have so far treated the scalar field perturbations 
generated during a de Sitter stage, when the evolution of 
the scale factor is a priori fixed independently of the 
behaviour of the scalar field. However, if this scalar 
field is the inflaton, it is the dominating component 
driving the accelerated expansion; its perturbations will 
then in turn modify the energy-momentum tensor, induc- 
ing fluctuations in the metric and specially in the scale 
factor evolution, via the slow-roll Friedmann equation 
H 2 ~ (8nG/3)V(4>). Besides scalar field perturbations 
Sip, we must now also consider the small metric perturba- 
tions. Both the metric fluctuations and the inflaton per- 
turbation Sip will assume different values depending on 
the choice of the coordinate frame. Therefore, in order to 
avoid this ambiguity and deal only with physical degrees 
of freedom, it is convenient to define the gauge-invariant 
comoving curvature perturbation 1Z = ip + H(5ip/4>) 
which measures the intrinsic spatial curvature on hyper- 
surfaces of constant time Q ■ Defining in conformal time 
z = a 2 (f>'/a' and u = —zlZ, the latter variable satisfies 
the equation of motion 



VV 



0. 



(V.l) 



Expanding the last term to first order in the slow-roll 
parameters e = —H/H 2 and rj = V" /3H 2 , one finds 



z 
z 



1 



(V.2) 



where v 



§ + 3e - rj. 

Therefore, in the slow-roll approximation the gauge- 
invariant normal modes itk satisfy the same equation of 
motion 111.9(1 . the only difference being in the definition 
of the parameter v labelling the solutions. We can thus 
apply also to u the stochastic formalism we developed in 
Sec. [n]for the case of a massive scalar field in a pure de 
Sitter background, obtaining 



(V.3) 



We can then apply to TZ the results derived for the power- 
spectrum of the perturbations of a test scalar field, con- 
cluding that also for the curvature perturbation we have 
^(fc) oc |/i(fc)| 2 . 




2.5 3 fer, 



Figure 1: Conditional power-spectrum obtained with the 
Gaussian window ljll.26^ in the limit r»/ri n — oo, for dif- 
ferent values of a varying from 0.1 (left curve) to 0.9 (right 
curve) . 



More precisely, recalling (fill . 9|> and taking the limit 
k\r\ <C <J < 1 (which is reasonably satisfied on cosmolog- 
ical scales) we get 



Vn(k) 



A 2 



-k\u\ J u (k\T.\)H u (k\T.\)WUk\T.\) 



iW a {k\rJ) (k\r\) 



\2tj-6e 



(V.4) 



If the boundary term in (|111.9|) can be neglected, the 
power-spectrum becomes 

V n (k)=A%W^k\n\)(k\r\)^. (V.5) 

thereby showing, since W 2 < 1, the presence of a blue- 
tilt on the largest observable scales. However, in the 
limit where a fc|r*| (since W CT (fc|T-^ | ) ~ 1) we recover 
the ordinary result 



V n {k) = A 



\2r;-6e 



(V.6) 



For generic values of a, we show in Fig. ^ the power- 
spectrum obtained with the Gaussian window function 
(|II.26|) . where the only approximation is that of consid- 
ering Ti n ^ 7"*- For this particular choice of the filter, we 
observe a blue tilt of the power-spectrum for k < 3<7a*H*. 

This blue tilt stems from the fact that a smooth win- 
dow function does not make a sharp separation in Fourier 
space but it gradually weighs the modes, allowing for 
a small low-frequency contribution to the fine-grained 
part of the field (in term of which the noise is defined) 
while depleting modes whose wavelength is immediately 
smaller than the cutoff scale. The colored noise origi- 
nated from such a window is thus able to generate fewer 
fluctuations than a white noise on scales slightly smaller 
than the comoving coarse-graining domain. 

As a consequence, since we impose the constraint that 
in our comoving patch of Universe the fluctuations can 
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Figure 2: . Angular power spectra of the CMB temperature anisotropies obtained with the Gaussian window 111.2611 and different 
values of the constraint time r». The cosmological parameters used in the computation of the radiation transfer functions are 
Q. b = 0.046, Q,cdm = 0.224, Q A = 0.730 and h = 0.72, corresponding to the best fit of the WMAP data 0. Dots represent 
the WMAP Binned Combined TT Power Spectrum, obtained from the LAMBDA website: http://lambda.gsfc.nasa.gov/ 



grow only after r», the scales that are leaving the hori- 
zon in the following few Hubble times receive less "ran- 
dom kicks" before freezing out than in the white-noise 
case. Therefore, the power-spectrum is a function of k 
smoothly interpolating between the values and 1 it as- 
sumes for small and large k, respectively. 



This power-spectrum can be used to calculate the 
CMB multipolcs predicted by a specific choice of the win- 
dow function W . Quite generally, we will expect to find a 
suppression of the lowest multipole, which are sensitive to 
a modification of the power-spectrum on this very large 
scales. However, in order to quantify this suppression one 
needs to choose the shape of the window function and 
the precise time t* at which the constraint is set. For the 
Gaussian window (|II.26(1 . the predicted CMB multipolcs 
are plotted in Fig.|21for different values of r*. The cosmo- 
logical parameter used to compute (with CMBfast [Tsl ]') 
the radiation transfer functions are the best fit of the 
WMAP data The significance of the low multipolcs 
suppression varies depending on the choice of the con- 
straint time: the effect is strong for t* ~ 3000/i -1 Mpc 
(corresponding to the present horizon size), while it be- 
comes fainter for earlier values of r„ and it is practically 
absent for r* ^ 9000/1" 1 Mpc. 



VI. CONCLUSIONS 

In this paper we calculated the power-spectrum of the 
comoving curvature perturbation 1Z in the framework of 
stochastic inflation with colored noise. We can conclude 
that a careful analysis of the stochastic behaviour of the 
fluctuations, for a physically plausible choice of the win- 
dow function W a (therefore excluding the step function) , 
leads to a curvature power-spectrum showing a relevant 
blue tilt for scales such that k ~ cra^H*, correspond- 
ing to physical lengths about a^ 1 times greater than the 
present Hubble radius. 

The possibility of finding a blue tilt on observable 
scales will then depend on the size of the coarse-graining 
domain at which thehomogeneity constraint is set. If 
aa*H„ is of the order of aoH (today's Hubble radius), 
then also the biggest scales inour observable Universe 
get tilted, while if it is much smaller this effect might 
be completely unobservable. We computed the pre- 
dicted angular power-spectrum C;'s of CMB tempera- 
ture anisotropies for a particular choice of the colored- 
noise window function, finding a suppression of the low- 
est multipoles whose importance depends on the pre- 
cise value of t*: this effect actally becomes more and 
more evident as t* approaches 3000/1" 1 Mpc (the size 
of the present Hubble radius), whike it is practically 
absent for t* ~ 9000/i -1 Mpc. This is an interesting 
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feature because the recent measurements of the CMB 
anisotropy performed by WMAP seem to give evidence 
of a lack of power on the largest observed scales with 
respect to the predictions of the standard inflationary 
scenario [ToT ]. In our approach, this deviation from the 
standard result of an almost scale-invariant perturbation 
spectrum is achieved without any new physical input in 
the theory, being merely a consequence of the choice of 
initial conditions once some over-simplifying approxima- 
tions (namely, the white-noise choice) are removed. We 
find it encouraging that no introduction of new physical 
ingredients is needed to explain the anomalous behaviour 
of the lower multipoles of CMB anisotropies. 
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